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It has been shown that superconducting domain walls in a model with £7(1) x 
Z2 symmetry can form long-lived loops called kinky vortons from random initial 
conditions in the broken field and a uniform charged background in (2+1) dimensions. 
In this paper we investigate a similar model with a hyper-cubic symmetry coupled to 
an unbroken £7(1) in which the domain walls can form junctions and hence a lattice. 
We call this model the charge-coupled cubic-anisotropy (CCCA) model. First, we 
present a detailed parametric study of the £7(1) x 7L<i model; features which we vary 
include the nature of the initial conditions and the coupling constants. This allows us 
to identify interesting parameters to vary in the more complicated, and hence more 
computationally intensive, CCCA models. In particular we find that the coefficient 
of the interaction term can be used to engineer three separate regimes: phase mixing, 
condensation and phase separation with the condensation regime corresponding to 
a single value of the coupling constant defined by the determinant of the quartic 
interaction terms being zero. We then identify the condensation regime in the CCCA 
model and show that, in this regime, the number of domain walls does not scale in 
the standard way if the initial conditions have a sufficiently high background charge. 
Instead of forming loops of domain wall, we find that, within the constraints of 
dynamic range, the network appears to be moving toward a glass-like configuration. 
We find that the results are independent of the dimension of the hyper-cube. 



I. INTRODUCTION 

One idea to explain the mysterious dark energy which appears to pervade the Universe 
is that it is due to a network of domain walls which has frozen into some kind of static 
configuration which is akin to a soap film (TJ |2]. Such a model predicts an equation of 
state with w = P/p = —2/3 and can be represented in cosmological perturbation theory 
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by an elastic medium with rigidity and a relativistic sound speed (see ref. [3] for a detailed 
discussion). An important question is whether such a static network can be created from 
random initial conditions. This would require the formation of either a regular lattice or a 
more amorphous glass-like state. 

In this paper we aim to develop further the idea that charge condensation onto the walls 
can be used to slow down and possibly even prevent the natural propensity of random wall 
configurations to collapse under tension, lose energy as fast as causality allows and scale in 
self-similar way with the number of walls being oc £ -1 . This might be possible in models 
where the domain wall forming fields are coupled to a complex scalar field whose potential 
has an unbroken £7(1) symmetry. 

A model with U(l) x Z 2 symmetry was discussed in ref. |3] in (2+1) dimensions. It 
was shown that long-lived loop configurations can be formed from initial conditions where 
the domain wall forming field is distributed randomly among the two vauca and the U(l) 
symmetric field is set to a uniform charge background. These loop configurations, known as 
kinky vortons [51 E], are stabilized by the Noether current and associated charge in a similar 
way to cosmic vortons [7]. 

The U(l) x Z 2 model does not allow for junctions between domain walls and one might 
think that they are necessary in order for the model to form a lattice or glass-like state. 
Junctions are present in the cubic anistropy model which is a model for an iV-component 
field whose potential has minima on vertices or the centre of the faces of a hyper-cube. The 
dynamics of this model have been be shown to lead to something close to the standard 
scaling regime within the narrow dynamical range possible in numerical simulations [8]. 

We will start by performing a detailed parametric study of the U(l) x Z 2 model. We 
will show that our previous results are largely independent of the way the initial conditions 
are set up and will investigate the scaling behavior as a function of the coupling constants. 
Importantly we find that the specific choice of the coefficient of the cross-coupling used in 
ref. [3] is a boundary between two regimes: phase mixing which takes place for lower values 
and phase separation for higher values. Only at the exact value used in ref. [3] does the 
phase condensation, necessary for the modified scaling behaviour observed, take place. We 
then go on to study a cubic anisotropy model coupled to a complex scalar field. We will call 
this model the charge-coupled cubic anisotropy (CCCA) model. Once we have identified the 
charge condensation regime, we find that this model behaves a similar way to the U(l) x Z 2 
model with charge becoming localized on the walls and we go on to discuss the nature of 
the configurations which are produced. 
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II. MODELS AND NUMERICAL IMPLEMENTATION 

The models we will investigate have two interacting fields $(a; M ) and a{x^) such that 
the equations of motion permit domain wall solutions in the $-field and a has a conserved 
Noether charge. We will consider models with a Lagrangian density which has the generic 
form 

c = l -d^d^ + l -d»ad,a - v, ($) - v 2 (H) - m VI 2 , (2.1) 

where the conserved charge of the complex scalar field a{x^) G C is given by 

Q = J d 2 xp Q = j i J d 2 x (aa-cra), (2.2) 
and Vi($) has discrete minima. We will take the specific form 

^(M) = ^(M 2 -^) 2 , (2.3) 

and the symmetry breaking term V]_($) to be one of two cases 

Vi a \$) = ^($ 2 -r/|) 2 , (2.4a) 

A N 
= -f (|$| 2 -4) 2 + ^^, * = (<j> 1 ,<f> 2 ,...,<f> N ),<l> i (af)eR. (2.4b) 



i=i 



We can rescale the energy and length units allowing us, without loss of generality, to define 
A$ = 77$ = 1. The model defined by is a special case of V} b) with e = 0, N = 1 and taking 
(3 = it becomes the non-interacting discrete Z 2 model which has a static kink (domain wall) 
solution. For (3 7^ 0, the model has superconducting domain wall solutions [5j 19] for specific 
values of the parameters. The model defined by V/ for e = has a global O(N) x U(l) 
symmetry. For (5 = the O(N) symmetry is broken to a discrete hyper-cubic symmetry 
group when e 7^ 0; for e < the points of the vacuum manifold are the 2N faces of an iV-dim 
hypercube, whereas for e > it is the vertices of the hypercube. For e > and N > 2 there 
are walls of a number of different tensions whereas only 2 different tensions exist for any iV 
with e < 0. 

The vacuum in the e > case in the presence of charge needs subtle consideration as we 
shall now describe. We shall continue with an iV = 2 case for simplicity, but the method can 
be generalized by the embedding Mr — » S 1 ^^ 1 . First let us write <f>\ = |$| cos#, 02 = \$\ sin# 
so that the N = 2 CCCA model's potential reads 

v = ^ (|$| 2 - vl) 2 + ^ (M 2 - v l) 2 + m 2 W\ 2 + e|$| 4 (cos 4 e + sin 4 e) . (2.5) 
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As all terms are positive definite (except the last term) and the (uncharged) vacuum has 

|$| 2 = -^L e>0, (2.6a) 

2e + A$ 

i*' 2 = « <»■ < 2 - 6b > 



the final term in (2.5) is only one which can select the discrete points that make up the 
vacuum manifold. It is now a trivial task to note that the vacuum for the e > case is 
constructed from the phase-set 9 > = {7r/4, 37r/4, 57r/4, 77r/4} and in the e < case by the 
phase-set 9 < = {0, ir/2, ir, 3tt/2} (which represent the vertices and faces respectively of a 
square). The height of a potential barrier between adjacent minima can then be computed, 
AV = V(6 e 0<) - V(6 e 9>), and is easily found to be AV = ^|$| 4 e. The strong 
dependence of the height of the potential barrier AV upon |$| means that if < 1 

only a very small perturbation is needed to surmount the potential barrier and move between 
minima. In the presence of coupling to charge (i.e. (3 ^ 0) this effect is so acute that regions 
of space occupying the discrete minima do not form in the e > case which has | $ | /r]& < 1 
in the vacuum (whereas |$|/^$ > 1 in the vacuum in the e < case) and hence in all results 
presented in this paper we use e < 0. 

In order for the discrete hyper-cubic symmetry to be broken locally one requires that the 
effective mass for $ is negative, m e 2 ff $ < 0. Including the interaction term corresponds to 
the constraint 

ml^ = P\°\ 2 - Y^l<0, (2.7) 
and hence one requires for symmetry breaking that \a\ 2 < k where 

K ^-f. (2.8) 

The fact that |cr| 2 varies from place to place means that the symmetry breaking condition 
m cS $ < is n °t necessarily met locally providing a rich landscape of phenomenologies. 

For our numerical work we evolve the equations of motion in Minkowski spacetime in 
(2+l)-dimensions by discretizing on a regular square grid of P grid-points along each side 
with grid-spacing Ax. We approximate spatial derivatives to fourth order and time evolution 
is achieved using a second-order leapfrog algorithm with time-step size At. We always use 
Ax = 0.5, At = 0.1 and we use both P = 1024 and P = 4096 - the former for computational 
brevity and the latter for high dynamical range. The boundary conditions we employ are 
periodic which sets an upper bound, r, on the time a given simulation can be used to deduce 
dynamical information. This is the time taken for an emitted signal to interact with itself 
after having crossed the periodic boundaries - the so called "light crossing time" . It is simple 
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to show that r = P Ax/2 At. Our results on the evolution of the scaling of domain walls are 
only valid upto this time; however, structures that persist after this time are of interest since 
they indicate stability of the wall network against perturbations due to radiation which is 
ubiquitous within the simulation. 

Typically we use initial conditions that could describe a system directly after a phase 
transition in a uniform charged background. To achieve this, we assign every grid-point a 
random position on the vacuum manifold (with $ = 0), and set a = Ae~ lujt . This corresponds 
to an initially homogeneous charge density Pq(0) = A 2 u specified by two parameters A, u. 
An unphysical side-effect of our initial conditions is that large energy gradients are setup 
between adjacent grid-points. To remedy this we introduce a damping term which acts on 
only the $-field for the first 200 time-steps in our simulation - this reduces the overall energy, 
aiding condensation into domains, whilst leaving the total charge constant. The dynamic 
range of a simulation is 200 < t < r: the time period over which physically meaningful 
results can be extracted. 

In order to investigate the scaling dynamics of the domain walls, we use an algorithm to 
compute the number of domain walls at each timestep. This works by deciding if adjacent 
grid-points occupy the same point on the vacuum manifold - if they do, then no wall exists, 
and if different the wall count is incremented. We compute scaling exponents by fitting a 
power-law to the number of domain walls, iVd w oc t~ 7 , in bins of 200 timesteps over the 
dynamic range. 



III. PARAMETRIC STUDY OF THE U(l) x Z 2 MODEL 
The aim of this section is to present results from a parametric investigation of the U(l) x 



Z 2 model using (2.4a) as the symmetry breaking potential. These results are an extension 
to those presented in our previous paper jl] allowing us to identify interesting effects in the 
computationally less expensive U(l) x Z 2 model which might be important in the CCCA 
model. Typically we will vary one parameter and keep the others fixed to fiducial values 
A = 1, u = 1, X a = 1, P = 1/2 and rj = I. 



A. Variation of initial conditions and r/ a 



We have performed a set of simulations varying the initial charge parameters A and u. 
Figure [T] shows the evolution of the number of walls as a function of varying initial charge 
density - where we have used A to vary Pq(0). These are equivalent to results presented in 
ref. jl]. As one can clearly see the evolution diverges from the standard scaling law as Pq(0) 
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Figure 1: Evolution of the number of walls for various Pq(0), where we have fixed to = 1, and 
varied A. Each plot is an ensemble average over 10 realisations with P = 1024. In all plots we use 
j3 = 1/2, Tfe = 3/4. Also shown for comparison is the oc t _1 scaling law adhered to by the 
Pq(0) = 0.0 system. 

increases; a conclusion which is further strengthened by considering the scaling exponents 
presented in Figure [2] Also shown in Figure [2] are scaling exponents for simulations where 
A and u are separately varied to give the same values of Pq(0). It is clear that, within the 
computed errors, the two methods give similar results; albeit those with varying to being 
marginally lower. 

In Figure [3] we present the evolution of the number of walls as we vary the parameter 77^ 
for Pq(0) = 0.25. Evidently, we see a deviation from scaling for the range of rf a presented. 
As r] 2 a is increased the scaling exponent 7 diverges from the scaling value 7 = 1. The known 
kinky vorton solutions exist only for rfc > 1/2 since the analytic solution on which they are 
based requires this constraint, and therefore it is interesting to see that it appears that a 
superconducting wall solution exists outside the analytically known range. 



B. Variation of f3 and phase separation 



The parameter whose variation produces the most interesting results is the cross-coupling 
constant /3 since by varying its value we can achieve three different kinds of behavior which we 
shall call phase mixing, condensation and separation. In order to understand this behaviour 
we have done three things. First we will present results from ID simulations which make 
it easy to see what is going on, some 2D simulations equivalent to those presented in the 
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Figure 2: Scaling exponents, where N^ w oc t~ 7 , for different values of /?q(0) created by varying 
either A or u), whilst keeping the other fixed. Each panel is for a specific value of pq(O), where the 
varied parameter is shown in the legend. All plots where u is varied are offset by 50 time-steps to 
allow the lcr-error bars to be seen. We have used P = 1024, rfc = 3/4 over 10 ensembles. 

previous section and finally some simple analytic arguments which suggest the existence of 
the phase separation regime. 

For our ID simulations the initial conditions have half of the grid with the $-field in one 
vacuum and the other half in the other vacuum and a homogeneous charged background. 
The state of the system is presented at t — 800 in Figure [4] for three different values of (3. 
For reasons that will become clear later, we define {3 = \\/\<$,\ a = \ and we will consider 
the cases (3 = ^(3 , (3 = (3 and (3 = |/3o. 

The case with (3 = |/3o has evolved to the standard kink solution and the charge density 
is homogeneously distributed. The case with (3 = (3$ has evolved to a superconducting kink 
solution with the charge density localized on the kink between adjacent vacuum states. Away 
from the kink the charge density is lower but non-zero which leads to the <3>-field being away 
from the vacuum points ±1. The final case, (3 = |/3 , leads to the production of two types of 
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Figure 3: Evolution of the number of walls and scaling exponents A^w oc t~ 7 for various rfe with 
P = 1024. Each case is an ensemble average over 10 realisations. In all plots we use A = 1/2, 
uj = 1, so that Pq(0) = 0.25 and /3 = 1/2. We have only presented a few of the scaling exponents 
of rfe, for clarity - the trend is clear. 
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Figure 4: Snapshot of the <j?-field and charge density pq at t = 800. The initial configuration has 
each half of the grid in each of the vacua with an initially homogeneous charge density with oj = 1, 
A = 1/2 (i.e. />q(0) = 0.25). As indicated, we display the results from (3 = ^/3o,/3o, |/5o; where 
/3o = \ ■ For /3 < (3q, one observes the homogeneous charge distribution and a domain wall only 
between vacuum states. For (3 = (3q, the charge distribution has evolved to localize on the domain 
wall separating vacuum states. For j3 > (3q, there is the production of a domain wall between the 
components of the mixture as well as walls between vacuum states; thus, domain walls interpolating 
between vacua $ = 1 and <3? = — 1, and between $ = and |$| = 1. 
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walls, one which is the superconducting wall and a second "domain wall" which interpolates 
between regions with and without charge - the so-called phase separation regime. We call 
the first regime (/3 < (3 ) the phase mixing regime, the second ((3 = (3 ) the condensate 
regime and the third {(3 > flo) the phase separation regime. Thus, we have a description of 
the system as being a mixture of two components: $ and a, and we can see that there are 
domain walls between states | ($) = 0, (er) ^ 0) and | ($) ^ 0, (a) =0), as well as between 
adjacent vacuum states |<& = +1) and |$ = —1). From here and on we shall call the wall 
between adjacent vacuum states a vacuum-wall and the wall between the components a 
component-wall. The vacuum-wall may also have a superconducting condensate localized 
upon it (as is the case in the /3 = (3q system) in which case it is a superconducting vacuum- 
wall. 

We have evolved the 2D equations of motion using (3 = |/3 , (3 , |/3 and we present 
results in Figure [5] We show the evolution of the $-field from random initial conditions 
and a homogeneous background charge. We have fixed = 3/4,pq(0) = 0.25. We see that 
the three regimes identified in ID realize themselves in 2D with the production of vacuum- 
walls and component-walls where appropriate. The evolution of the (3 = \(3q case proceeds 
unimpeded by the charge which is confirmed by the computed scaling exponent 7 « 1. As 
is increased towards (3 , the scaling exponent 7 — > 0. For the case of j5 = /9 there are 
vacuum-walls which clearly do not collapse in the standard way - upon inspection of the 
other fields, we find high charge and current densities along the domain walls as explained 
in ref. jl]. We have argued that the superconducting condensate has become naturally 
attracted to the vacuum-walls and the conserved current provides a restorative force which 
opposes the natural tension in the vacuum-walls. 

The third case, (3 = |/5o, shows a very different picture. One can clearly see that regions 
where $ occupies the discrete vacuum states (i.e. where $ = ±1) appear from a background 
of |$| ~ 0. Upon inspection of the charge density, we find that |pq| « inside a region 
occupying a vacuum state and |pq| peaks in the regions where |$| ~ 0. These observations 
corroborate those of the simple ID system (Figure [1]) and confirm the veracity of our classi- 
fications of phase mixing and separation. We rely on the visual observation that the domain 
wall network has not entered the standard scaling dynamics in the phase separation regime 
since the algorithm for calculation of N^ w is no longer valid. If we let \(3 — (3o\ — > in the 
phase separation regime (i.e. approach from above) then separation of the components still 
occurs but the extent of separation decreases as \(3 — (3 \ — > 0. 

One can understand the behavior observed in the simulations using simple analytic argu- 
ments. We will show that a component-wall exists between the $ and a fields, for (3 > 
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Figure 5: Evolution of the $-field for j3 = hftoi Po, |/?o (top to bottom) at times t = 
400,800,1600,3200,6400 (left to right) with p Q (0) = 0.25,7/2 = 3/4 and P = 1024. Each 
colour /shade represents each of the two points on the vacuum manifold (with purple/darkest 
being |<3?| ~ 0). The behavior observed in these three cases is similar to that presented in ID 
in Figure [4j Vacuum- walls are observed when regions in different vacuum states are adjacent and 
component-walls where regions of vacuum occupation are adjacent to regions of extended charge 
density. 

in the U(l) x Z2 model. In particular, we substitute the ansatz for a y-directed vorton 

(Q\ ( cos(£/2) \ 

\<t) P \sm(6/2)e i ^+ k y)J , { } 

into the energy density where 9 = 9(x) and p = p(x) is the total field density p 2 = $ 2 + \a\ 2 . 
Assuming A$ = A CT = A we find that we can write the free energy 

F = E-pQ = |d 2 x{^ 2 |VC + Ap 4 sin^ + ^|Vp| 2 + ^(p 2 -^ + ^4| ; (32) 

where the charge and chemical potential are denned to be 

Q = iu J d 2 xp 2 sm 2 (6/2), p = ±- [A (tjI - r/ 2 ) + u 2 + k 2 ] , (3.3) 

and A = i(/3-iA). 



11 



For A ^> |A|, p will be approximately constant and hence 

F^Jd 2 x j^p 2 |W| 2 + Ap 4 sin 2 #J. (3.4) 

This is the double sine-Gordon model and is known to have solutions for sufficiently high Q 
if A > [10]. The solutions have 9{x) = tt for \x\ <C R and 8(x) = for |x| ^> R, which 
is a region of charge with a domain wall separating the two phases. We note that A = 
corresponds to (3 = (3 = |A and hence this phase separation regime is valid for (3 > (3 Q . 



C. Zero Net Charge 



All the simulations we have presented so far have an overall net charge within the simula- 
tion box. Although it might be possible for some overall net charge to exist in the Universe, 
a more conservative idea would be that there is a zero net charge with regions of positive 
and negative charge having some correlation length different to that of the vacuum-wall 
forming field. In order to see if the phenomena observed in the case of a homogeneous 
charge persist in this case we have completed a set of simulations with P = 4096 and an 
initial charge-correlation length which is less than the box size but greater than that for the 
vacuum- wall forming field (which is just the grid spacing). We will set /3 = /3 and the net 
charge Q = but with a local charge density pq ^ 0. The initial conditions have alter- 
nating square regions of charge, each region being 1024 grid-squares across - every region 
has the same magnitude of charge density and adjacent regions being opposite in value. 
Thus, along a particular horizontal or vertical line, this amounts to the charge density being 

H 1 — . Images of the charge density in these simulations can be found in Figure [6] for two 

initial maximum charge densities. Naively, one may think that charge will find anti-charge 
and annihilate. Initially, some annihilation does occur on the boundaries between regions 
with different charge. However, we find that once charge has condensed onto vacuum-walls 
it stays there to create a superconducting vacuum-wall. One observes that vacuum-walls 
have both positive and negative charge condensed on them, as the vacuum-walls pass be- 
tween regions of different charge densities. The resulting superconducting vacuum-walls do 
not collapse as readily as their uncharged counterparts - although the scaling exponent 7 
is not as far from unity in this correlated-charge case as in the homogeneous charge case: 
(1 — 7) cor < (1 — 7) homo- This is due to some charge density annihilating reducing the 
amount of charge available to condense onto the walls. 



12 




Figure 6: Evolution of the charge density pq in simulations of the U(l) x Z2 model in the condensate 
regime. Initial conditions have charge density /Oq(0) = 0.09,0.25 (top to bottom) in blocks as 
described in the main text. Images are at t = 3200,6400,12800 (left to right) and are coloured 
such that red/blue denotes positive/negative and grey where the value is less than 10% of the 
maximum. We have used P = 4096 and rfe = 3/4. 

IV. GENERALIZED PHASE MIXING AND SEPARATION 

In the i7(l) X Z2 model the parameter /3 = 0o, the value of the cross coupling above 
which one is in the phase separation regime, corresponds to the determinant of the quartic 
coupling matrix being zero. The concept of this "mixture stability" occurs in a number of 
physical systems. The simplest example is that of oil and water where bubbles of oil separate 
out from the water and are stable to decay essentially due to the system being in the phase 
separation regime. It has also been observed in two component Bose-Einstein condensates, 
for example, between the hyperfine states of Rubidium-87 (TTJ H2], and can be useful in the 
creation of topological defects [13J. In this section we will derive the equivalent condition 
for the CCCA model which is a little more complicated due to the additional fourth order 
terms proportional to e. This will allow us to identify the value of /3q in these models and 
test the concept in subsequent sections. 
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We derive the mixture stability criterion by computing conditions on the coefficients such 
that the stable configuration can have both components in mutual existence. To begin, let 



us write the potential term using (2.4b) as 
1 1 



N 



V = -£n|$| 4 + -£2 2 M 4 + 0i 2 |$lV| 2 - mi|$| 2 - m 2 \a\ 2 + e^0 4 + V , 



(4.1) 



i=i 



where we have introduced interaction coefficients gij and masses m*; this is to make contact 
with the notation used in the condensed matter literature. The gij,rrii and Vq are given in 
terms of the model parameters by 



9ll — 2^*' 922 — 9l2 — Pi 

m 1 = \X^r]%, m 2 = ^rfc, V = \ (A$ry| + Kvt) ■ 
In addition, we introduce the number density parameters, 



n<s> 



N 

E 

i=l 



Tl4 « 71/4 



n„ 



(4.2a) 
(4.2b) 



(4.3) 



so that the potential (4.1) is given by 



1 1 N 

v = ~9nn% + r</22"ff + 9\2n^n a - m x n^ - m 2 n a + e n? + V c 



(4.4) 



i=i 



To find the field configuration that extremizes the potential, one computes the first deriva- 
tives of the potential with respect to each of the density parameters (each ^-component to 
be computed separately) setting them to zero. This implies that 

dV 
drii 
dV 



g u n^ + g 12 n a - m 1 + 2erii = 0, 
022?v + g\ 2 nq, - m 2 = 0. 



(4.5) 
(4.6) 



As each of the N components of (4.5 ) are individually zero, so is their sum. Hence, summing 
(4.5 ) over i — 1, . . . , N gives 



dV 



E 



dV 

drii 



(<7nn$ + 9\ 2 n a - m^N + 2en$ = 0. 



(4.7) 



This step of constructing an equation in terms of n$ and n a without reference to the individ- 
ual Hi is crucial in what follows - not all models can be written in such a way, for example, 
if the anisotropy term has non-diagonal quartic terms. 

The condition that a given configuration is stable comes from the requirement that the 
eigenvalues of the Hessian are positive-definite. That is, solving 

d 2 V 



dn a dn/3 



A5, 



a/3 



with n a G {n$, n a } } 
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and requiring that A > 0. This results in the condition 

fan + 2e)g 22 - g\ 2 > 0. (4.9) 
If e = the left hand side of the above inequality is simply the determinant of the matrix of 



quartic coupling coefficients between the components of the mixture. Finally, (4.9) in terms 



of the model parameters (4.2) is given by 



P<Po, f3 = \V^K + ^K- (4-10) 

This is the condition that must be satisfied for a composite mixture of $ and a to be stable. 
Hence, for > 0q, the mixture is unstable, with the components separating. Figure [7] has a 
schematic describing the regimes discussed here. 

condensate regime 
phase mixing phase separation 



/3 < A) P = A) p> A) 

Figure 7: Schematic depicting the parameter ranges for the three phase-interaction regimes of the 
17(1) x Z 2 model. 

V. CHARGE-COUPLED CUBIC ANISOTROPY (CCCA) MODEL 

Let us now turn our attention to the CCCA model. There is much the same phenomenol- 
ogy as in the U(l) x Z 2 model, the main difference being the presence of multiple vacua, and 
hence the possibility of junctions for N > 1. Throughout we take A CT = 1, r\ a — \/3/2 and 
e = —0.1. As with our U(l) x Z 2 simulations we present evolution of a random distribution 
of vacua in the <£>-field in a homogeneous charged background for a variety of values of the 
coupling parameter, 0, and initial charge density Pq(0); we vary A and fix uj = 1. 

To begin, we vary either side of 0q as we did for the U (1) x Z 2 model and take N = 2. 
Figure [8] shows the evolution of the (</>i, 2 )-fields where (3 = ^/3 , fa, |/3 . The case of 
/3 = ^/?o is very similar to that in the U(l) x Z 2 model: vacuum- walls evolve with scaling 
exponent 7 « 1 and the background of charge remains approximately homogeneous which 
is indicative of the phase-mixing regime. The condensate regime = /3q clearly does not 
evolve as quickly, with charge and current becoming associated with the vacuum-walls which 
appears to inhibit their collapse. It is interesting to note the existence of both X- and V-type 
junctions, especially at late time (i.e. t = 12800). The final case = |/3 is qualitatively 
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Figure 8: Evolution of the (0i, ^-fields in the regimes (3 = ^ /?o , /3o , § /3o (top to bottom), at 
times t = 800,1600,3200,6400,12800 (left to right), where p Q (0) = 0.25 and P = 4096. Each 
colour /shade represents one of the four points on the vacuum manifold with black/darkest being 
|<J>| ~ which corresponds to regions of high charge density. 

very similar to that in the U(l) x Z 2 model: it appears that the system is in the phase 
separation regime and we observe component- walls between regions with charge and those 
without. 

Figure [9] shows the configurations of the potential energy density ppe, real part of the 
condensate field 9?(cr), charge pq and current |J| 2 at t — 6400 of the (3 = f3 simulation. 
Figure [9|a) is not surprising: potential energy peaks on vacuum- walls as the field interpolates 
over maxima of the potential. Figure ^b) shows that 3?(cr) varies from positive to negative 
along the walls which indicates that the cr-field winds along a vacuum-wall - the alternating 
red-blue along a given vacuum-wall. From Figure and (d) it is obvious that both 
charge and current have become almost exclusively associated with the vacuum-walls. The 



inspection of an equivalent set of simulations presented in Figure [TO] for a lower initial 
charge density Pq(0) = 0.09 show more clearly the winding. It is the association of charge 
and current with vacuum-walls which provides the stabilization mechanism of kinky vortons 



[I] and it appears from Figure M and Figure 10 that we can extend the mechanism to the 



CCCA model. If we increase the threshold used to create the images of the charge density 
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Pq (so that the charge density at a given location must be a larger fraction of the highest 
value than that used in Figure [9] and Figure 10) it is revealed that junctions carry a higher 
charge density than vacuum- walls, PQjunc > PQ.waii- 

The images in Figure 11 show how the pseudo-stable structures evolve for times beyond 
the light-crossing time t ^> r in the regimes (3 = (3 , |/3 (note that P = 1024 was used 
to reduce computational time). As time progresses in the case of (3 = |/3o the separation 
between the vacuum-occupying and charge components becomes more obvious. At late times 
the configuration with (5 = (3q appears to have stabilized into a glass-like state. 

We have performed a number of simulations varying the initial charge density Pq(0) and 
the cross coupling parameter 0. The results are in Figure 12 at t = 3200 for P = 1024. 
These clarify how the system responds to different amounts of initial charge in the phase 
separation regime. Interestingly we find that, in the (3 = \(5q simulations, all images are 
qualitatively identical with 7 w 1 irrespective of Pq(0). Finally, as the initial charge density 
is increased in the condensate regime, it is obvious from the images that the vaccum- walls 
retain their length for longer. 

We have also performed a number of simulations for different values of N with /5 = /3 and 



Pq(0) = 0.25. Evolution of the fields are presented in Figure 13 and plots of the evolution 
of low and high tension vacuum- walls in Figure 14 It appears that the evolution for higher 
values of N is qualitatively the same as for N = 2, although the number of high tension 
walls decays quicker in higher- N models. This is probably due to the larger phase-space of 
low tension wall trajectories that high tension walls can decay into for large N models. It 
is clear from the plot of the evolution of the number of walls that the systems have scaling 
exponents 7 < 1. 



VI. CONCLUSIONS 

In this paper we have developed the concept of charge condensation onto a random 
domain wall network and we have given analytic and numerical arguments for its causes 
and consequences. Furthermore, we have shown that random domain wall systems can 
form a pseudo-stable glass-like network with the key addition to previous models being the 
interaction of the field with discrete minima and a field with conserved continuous U(l) 
symmetry. Because conserved symmetry currents condense on domain walls in a particular 
parameter regime (what we called the condensate regime) we can extend the mechanism 
that stabilizes kinky vortons to the rather more general CCCA model. 

The justification for this comes from the observation that in the phase mixing regime 
the charged field remains homogeneous and does not condense on domain walls, which 
promptly collapse under their own tension and the network enters a standard scaling dy- 
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Figure 9: Properties of the field configurations at t = 6400 for /3 = /3q with N = 2,P = 
4096, pq(0) = 0.25. The colouring scheme is such that red/blue denotes positive/negative and 
grey being 9?(cr) < 40% maximal value; /Ope,q < 10% maximum; |J| 2 < 2% maximum. These 
images show that charge and current have become almost exclusively associated with the vacuum- 
walls. 
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Figure 10: The same as presented in Figure [9] but for a lower initial charge density, /Oq(0) = 0.09. 
The formation of junctions and the fact that the <r- field winds around the walls is more obvious at 
this lower charge density. 
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Figure 11: Evolution of the (<^>i, $>2)-fields in the regimes j3 = /3q, |/3o (top to bottom), at t = 
t, 2t, 3t, 4t (left to right), where pq(0) = 0.25 and P = 1024. One can observe that structures 
have persisted beyond the light crossing time r. In the /3 = /?o regime the network's stability 
is due to symmetry currents generating a superconducting vacuum-wall solution which freezes 
the network into a glass-like lattice. In the j3 = |/?o regime phase separation has prevented two 
adjacent vacuum occupying regions meeting and annihilating (as is the case in the phase mixing 
or uncharged systems). 



namic with A^ w oc t~ 7 ,7 = 1. However, in the phase separation regime the formation of 
a component- wall appears to prevent vacuum-walls from touching each other and annihi- 
lating, thus preventing the natural collapse of wall networks. In the condensate regime we 
have shown evidence for a winding of the complex scalar field along vacuum-walls and the 
association of the conserved charge with the wall network. We have shown that the phase 
separation phenomenon exists in both the U(l) x Z 2 model and the CCCA model. In both 
the condensate and phase separation regimes a non-standard scaling behaviour is observed 
with 7 < 1 and the extent to which a network stabilizes is controlled by the value of the 
initial charge density. In addition we have also shown that there does not need to be a 
net charge in a simulation for this stabilization mechanism to occur, merely a local charge 
density with a sufficiently large correlation length 

Controlling whether or not a theory is in the phase separation, phase mixing or condensate 
regime is achieved by the value of the cross coupling parameter (3, where a given theory has 
a fixed value (3q given in terms of the other quartic coupling parameters. When (3 > (3 the 
theory is in the phase separation regime, when /3 < (3 Q the theory is in the phase mixing 
regime and when /3 = flo the theory is in the condensate regime. 
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pq(0) = 0.09 p Q (0) = 0.16 p Q (0) = 0.2 p Q (0) = 0.25 p Q (0) = 0.36 



Figure 12: Snap-shots from the evolution of (<^i,<fo) at i = 3200, with P = 1024, in (/5q(0),/3)- 
space. Note that only vacuum-walls exist for /3 < /3q and the existence of component walls for 
f3 > /3q. Regardless of the value of Pq(0) in the regime (3 = ^A) the network has collapsing 
vacuum- walls, with scaling exponent 7 ~ 1. One can observe that as pq(0) increases in the regime 
j3 = fto, the glass-like nature of the resulting lattice becomes more apparent. As pq(0) is increased 
in the /3 = |/3o regime, there is more charge to separate from the vacuum-occupying regions, which 
is obvious by noting that the size of the charge-regions increases with initial charge density. 

The dynamics of wall networks was studied extensively in refs. [1~4T[TT] . On the basis of 
these works the authors suggested a "no frustration conjecture". They claimed to prove this 
based on their so called "ideal model" which has junctions and walls of identical tensions. 
The studied scaling regime occurs because the network loses energy as fast as causality will 
allow; there is nothing to prevent this. The results presented here suggest an important 
caveat to this: the no frustration conjecture will be true only in the absence of charge, or 
some other phenomena which prevents the collapse of walls. 

We have shown that it may be possible to stabilize networks. Hence, we have provided 
theoretical evidence which might lead to a physical manifestation of the elastic dark energy 
model, which, until the charge-coupling mechanism presented here, had been hampered by 
the oc t~ x scaling behavior of domain wall networks. However, still more works needs to be 
done before we can claim to have a fully consistent framework for domain walls to provide 
the dark energy of the Universe. 
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Figure 13: Evolution of the ^-fields for N = 5,4,3,2 (top to bottom), at t = 
400,800,1600,3200,6400 (left to right), where p Q (0) = 0.25, /3 = (3 and P = 1024. As the 
number of vacua is increased with N, the network dynamic is qualitatively identical to the N = 2 
case where charge and current form superconducting vacuum-wall solutions stabilizing the network 
into a glass-like lattice. 
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